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Abstract 

We study the asymptotic behavior of a diffusion process with small diffusion in 
a domain D. This process is reflected at dD with respect to a co- normal direction 
pointing inside D. Our asymptotic result is used to study the long time behavior of 
the solution of the corresponding parabolic PDE with Neumann boundary condition. 
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1 Introduction 



Consider the following parabolic initial-boundary value problem 
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u%x,0) = g{x) , 



e>0 ; 

x(^DV^^D ; 
xedD , t>0 . 



(1.1) 



Here D is a d-dimensional bounded domain in M with smooth boundary dD. The 
initial condition g{») is smooth in Z? U dD. The matrix a{x) = (aij(x))i<jj<rf is posi- 
tive definite. The functions aij{x) are smooth and uniformly bounded, with uniformly 
bounded derivatives. There is a constant 9 > such that for any ^ = {^^,...,C ) we 

d 

have 02|^|2^ < J2 a^j{x)C^^ < O'^lClla- The vector field 6(x) = {b\x), ...,b'^{x)) have 

terms which are uniformly bounded, smooth in [D] (here and below [D] is the closure 
of D in Euclidean metric), and have uniformly bounded derivatives. The vector field 
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7(x) = (7^(2;), ...,7'^ (x)) is the inward co-normal unit vector field on dD with respect 
to the matrix a~^(x) = (a*-'(x))i<jj<rf = {aij{x))^^- ■<^. That is to say, for any vector 

d .. 
v{x) = {v^{x), ...,v'^{x)) tangent to dD we have (7, f)a-i(^) = J2 a^-' {x)'^^ {x)v^ {x) = 0. 

(Here and below (7, v)a-i(a;) is the inner product with respect to the matrix a^'^ix). For 
a detailed discussion of the co- normal condition we refer to [31 Section 2.5].) We also 
have |7|Rd = 1. 

Let us assume that the vector field h{x) is pointing outward to D on a connected 
subset diD of dD, and it is pointing inward on d2D = dD\diD. (It is never tangent 
to dD.) Let b{x) be the field coinciding with b{x) everywhere except at those points 
of diD. At these points b{x) is defined as the projection of b{x) onto the direction of 
the boundary. Suppose the dynamical system xt = b{xt) has all its w-limit sets on diD. 
These w-limit sets are points 0\, ...,Oi {I > 1). 

Our goal in this paper is to describe the long-time behavior of the solution u'^{x, t) 
of (1.1) as e —7- and t —> 00. One can relate problem (1.1) with a certain diffusion 
process X^ with small diffusion and refiection with respect to 7 on dD. This process 
can be described as a solution of the following stochastic differential equation: 



dXf = b{Xl)dt + ea{Xt)dWt + laD{Xl)^{Xt)det , Xq^ = x , CE = . (1.2) 

Here (t{x) is a (i x d matrix with smooth terms (and bounded derivatives) that 
satisfies a{x)a'^{x) = a'^{x)a{x) = a{x). The function 1q£){») is the indicator function 
of dD. The processes Xf and ^f are continuous time stochastic processes, adapted to 
the filtration {J-'t)t>o- They satisfy the following assumptions with probability 1: 

(1) The process Xf e [D] ; 

(2) The process ^f is non-decreasing in t and increases only at A = {f ; X^ € dD}; 

(3) The set A has Lebesgue measure zero. 

Under these assumptions, it was proved in [I] (also see [TTj) that such a pair of 
processes (Xf,^f) exist and is unique (in the sense of probability 1). The process ^f is 
called the local time of the process Xf on dD. (We remark here that this notion of the 
local time for the multidimensional diffusion process extends the classical 1-dimensional 
local time in [8j. See [IT] for a discussion based on SDE approach. For other discussions 
of the local time for multidimensional diffusion process we also refer to [9] and [IQj.) 
The process Xf is a strong Markov process in [D] and it satisfies the Doeblin condition, 
which leads to the existence and uniqueness of an invariant measure in [D] . 

It turns out that the solution u'^{x,t) of (1.1) can be represented as u^{x,t) = 
Eixg{Xf;) (see Section 4 for details). Thus the asymptotic behavior of solution u^{x,t) 



as e — > 0, t — > cx) is determined by the asymptotic behavior of the process Xf. However, 
the latter can be calculated using the Freidlin-Wentzell large deviation theory (see [5], 

ID)- 

In Section 2 of the present paper we will give an expression of the action functional 
S^rp of the process Xf . By using the large deviation principle for the family of processes 
{X|}e>o we will give a description of the asymptotic behavior of Xf in Section 3. Since 
the proof is based on the method of [SJ Ch.6] and [1], we will only prove some key 
technical lemmas and sketch the result. In particular, we give the algorithm on the 
calculation of metastable states. Section 4 provides the corresponding result for problem 
(1.1). We point out that a related question for elliptic boundary value problems was 
already considered in [B] (also see O Section 10.3]). An example is given in Section 5. 

2 Calculation of the action functional 

In this section we give an expression of the action functional corresponding to the 
large deviation principle of the process Xf. The main proofs and justifications of our 
results are contained in [1] (also see O Section 10.3]), so we just summarize the results 
we need. 

In [T], the authors have constructed the process {Xf,Q) corresponding to (1.2) by 
first realize it in the space M^ using the following stochastic differential equation: 

dy/ = 6(r(y/))(it + e(7{T{Yt^))dWt , Y^ = x^n% . (2.1) 

Here F : C[o,oo)(l^'^) ~^ ^[Q,oa)(^+) is a functional defined by 

r(V't) = (r(^))t ^ Tt{^) = {^l - A inf >1,^?, ..., V-f) (2.2) 

\)<s<t 

for ipt = {i^ti ■■■^i^t) ^ C[o^oo)(J^'^)- It was proved in [1] that in the case of a half space 
M^ one can take (X^^f) = (r(y/), (r(y/) - Y^Y). 

In the general case when D is a bounded region in W^ with smooth boundary one 
can take a finite covering of Z) by a set of open neighborhoods \Ui^...,Un}- Within 
each Ui {i = 1, ..., N) the process can be constructed via a homeomorphism between lAi 
and M , or between lAi fl D and M'J_ (when lAi fl dD 7^ 0). In the latter case we use the 
construction of the process in half space as above. By appropriately "glue" these pieces 
of the trajectories together one can construct the processes (Xf,,^f). The process Xf is 
the diffusion process with reflection in D and the process ^f is the local time on dD. 
For details of this construction we refer to [I], O Section 1.6]. 



It was shown in [H Section 1.2] that the corresponding action functional for the 
family of processes {Xf }e>o as e J, is given by the formula 



S±iip) 



1 f^ 

■^° (2 3) 

for if G Cjo^riCi-D]) absolutely continuous , ipo = x ; ^ ' 

+00 , for the rest of (/? G C[o^t]([-^]) • 



Here 

^(S = II I ^||2 V , 

and ||t;||a-i(a:) = (^'^)a-i(x) ^°^ vector V G M"^. 
We have 



= \qbWs)i{'^s) II . .||2 II , ^||2 V 

= -\qd{}Ps)t—, — r]j| [0A(6(c^s),7((^5))„-i(^^)] for a.s. sG [0, T] . 



Define 



h[x) = h{x)-\dD{x)^—^ [0A(6(x),7(x)),-i(,)] . 



(2.4) 

We see that h{x) is the field coinciding with h{x) everywhere except at those points 
of d\D. (Recall that d\D is the part of the boundary dD on which h(x) is pointing 
outward.) At these points h(x) is defined as the projection of h(x) onto the direction of 
the boundary. The action functional for the family of processes {Xf }£>o can now be 
formulated as 

5-+ ( \ ^ ] 2 / "'^'^ ~ K'Ps)fa~^{^,)'^^ ' for v? G C[o_r]([-D]) absolutely continuous , (/?o = a; ; 

[_ +00 , for the rest of (^ G C[or]([-D]) • 

(2.5) 

The deterministic trajectory Xf at which the above action functional is is also 

calculated in [Ij. It is given by the system xt = h{xi), xq = x, i.e., it coincides with the 

deterministic trajectory given by the vector field h{x) everywhere except at those points 

of d\D^ and at points of d\D it follows the projection of h{x) onto the direction of the 

boundary. 



We formulate below the large deviation principle for the family of processes {Xf }£>o. 

Theorem 2.1. (Large deviation principle) For the process Xf , we have 
(i) The set ^{s) = {ip £ CtQ^rp^dD]) : S'^((^) < s} is compact for every s > 0; 
(ii) Given (p G Cro7-i([D]). For any 5 > and any 7 > there exist an Eq > G such 
that for any < e < eo we have 

P{pot(^", v)<5}> exp[-e-2(5+,(^) + ^)] , (2.6) 

where T > and pori*, •) denotes the uniform distance between functions in Crg^Ti ([!?]); 
(Hi) For any 5, 7 > and any s > there exists an Eq > such that for any 
< e < £0 we have 

P{pot{X', ^s)) >5}< exp[-e-2(s - 7)] , (2.7) 

where Pot {(p,^{s)) = inf poTi^,ip)- 

3 Asymptotic behavior of X^ 

3.1 Estimates on the time to converge to oj-hmit sets on the boundary 

We now begin our study of the asymptotic behavior of the process X^. First, since 
the dynamical system xt = b{xt) does not have any oj - limit set within D, we shall 
expect that as e is small, the trajectories of Xf come to the boundary diD within finite 
time. (Notice that at points of diD the vector field b{x) is pointing outward and at 
points of d2D it is pointing inward. Therefore the deterministic trajectory X^ will not 
come to 820.) 

For any x,y £ [D], we define 

V^{x, y) = inf {S'Jy((^), ipo = x,ipT = y,^t ^ F> U dD, < t <T < 00} . 

VGC[u_t]([D]) 

Recall that the dynamical system xt = b{xt) has all its w-limit sets on diD. These 
w-limit sets are points Oi,...,Oi (/ > 1). Let us suppose, that for any x and y in [D], 
X ^ y we have at least one of V^{x,y) and V^{y,x) being > 0. 

For each Oi,i = 1,2, ...,/, by an a- neighborhood £a{Oi) of Oj, we refer to the 
intersection of D with an open ball having center Oi and radius a > 0. We use the 
symbol d£a{Oi) to mean the intersection of [D] with the boundary of the open a-ball 
centered at Oi. We call d£a{Oi) the boundary of the a- neighborhood of Oj. Let us 
choose a > such that the a-neighborhoods £a{Oi) for all Oi,i = 1,2,..., I, does not 
intersect each other. We now prove the following: 



Theorem 3.1. There exist positive constants c and Tq such that for all sufficiently 
small e > and any x € [D], Xq = x we have 

PxiCa >T}< exp[-e-2c(r - Tq)] , 
where (a = mf{t : Xf G U [^a(O^)]} • 

Proof. We consider the dynamical system xt = b{xt) on the whole domain DUdD, 

where vector field b{x) is defined as before. Since system xt = b{xt) does not have any 

oj - limit set in D and {b{x),^{x))a~'L(x)\d2D > 0, we can say that the time Ti(x) that 

the trajectory xt{x) spends until reaching diD is finite (if x £ diD, let Ti(x) = 0). 

Let y{x) be the point where trajectory first hits diD. Starting from y{x), the time 

T2{y{x),a) = T2{x,a) that the trajectory of system xt = b{xt) on diD spend to come 

I I 

into y [(£'2l(0j)] is also finite (as is the same, if y{x) € |J [£^2L(Oi)], then T2{x,a) = 0). 

The function T(x, a 
we have lim. 



Ti{x)+T2{x, a) is upper semi-continuous in x (i.e., for x, xq € [D] 
.x^xoT{x,a) < T{xQ,a)) because xt{x) depends continuously on x. Thus 

there exists Tq = max T{x, a) < oo. The set of functions in C ro,Tol ( [-^1 ) assuming their 

xe[D] 

values in [D]\ I |J £a(^Oi) ) is closed and thus Stj, attains a minimum A on this set. 



\i=l 



Taking into account the construction of Tq and the form of S^j, in (2.5), we see that 
A > 0. Let < (5 < f. Let $^(^/2) = {99 € C[o,To](^), V'o = a;, 5+,,, (99) < A/2}. We 
see that trajectories for which C,^ > Tq are at a distance > 6 from ^xi^f^)- Thus by 
the part (iii) of the large deviation principle we have 

PxiCa > To} < exp[-e-2(A/2 - 7)] 

for some < 7 < A/2. 

Thus by strong Markov property, 

PxiCa >{n + l)To} = I^xiCa > nTo; Px-^ {Ca > To}] 

< PxKa > nTo}exp[-e-^A/2 - 7)] . 
So by induction we see that 



Px{Ca>T} <P^{Ca> 



T 



To} 



<exp{-e~2(A/2-7) 



T 

Th 



<exp{-e-^(--l)(^/2-7)}. 



A/2 — 7 
Putting c = — — , we get as desired. D 



3.2 Transition probabilities between neighborhoods of the Oj's 

In this section we study the asymptotic transition probabilities between neighbor- 
hoods of the w-hmit sets {Oi, ..., O;}. We first provide several auxiliary lemmas. 

Lemma 3.1. There exists a constant L > such that for any x,y £ [D] sufficiently 
close to each other, there exists a function if G Cro^T]([-^])) V'o = x,(pT = y, such that 
ot( 



we have S^rp^ip) < L ■ \x — y 



Proof. Let x and y be so close to each other that they can be covered by one 
coordinate chart U. Let this coordinate chart correspond to a coordinate function 
li : [/ — 7- M ( or R'J_). The function u is smooth with bounded derivatives. Let us 
take T = \x — y\^d, 

ift = u'^ {u{x) + —{u{y) - u{x)) 
We have, for some constant M > 0, 



1 '■^ 






2 JO 



2 rT 



9 f 

< y / I'Ps - b{ips)\l,dds 



'0 
32 f-T 



2 







-{u ^y iu{x) + -{u{y) -u{x))j o{u{y) - u{x)) - b{ips) 



ds 



<e^M\hy-x\l,+Tma^\b{x)\l 
Taking into account that T = \y — x\^d, we are done. D 

Lemma 3.2. For any 7 > and any compact subset K C [D] there exists Tq such 
that for any x,y G K there exists a function '~pt,0<t<T,ipQ = x,ipT=y,T<TQ such 
thatS^r^iip) <V+{x,y)+j. 

Proof. We choose a finite (5- net {xi} of points in K; we connect them with curves 
at which the action functional assumes values differing from the infimum by less than 
^ and complete them with end sections using Lemma 3.1: from x to a point Xj near x 
and then from Xj to a point Xj near y, and from Xj to y. By choosing 6 small enough 
we get as desired. D 

We define 



F+(0„ O,) = inf {S+^{^) : ifo = Oi, ifr = Oj,^t (^[D]\ M {OJ, < t < T} . 

v'GCfo.T] ([£>]) jr. 

A "y+(Oi, Oj) version" of the above Lemma can proved similarly: one can take the 
curve if in such a way that it avoids Us^^i ji^s} and such that S'ot(v') ^ ^^(Oj, Oj) + 7. 
We omit the proof. 

Let constant po > be small. Let constant < pi < po- We denote by C the set 
D U dD from which we delete the po-neighborhoods of the Oi,i = 1,2, ...,/; by Fj the 
boundaries of the po-neighborhoods of Oj: Fj = d£pf^{Oi); by ^(j the pi-neighborhoods 
of the Oi , and by g the union of all the gi . 

We introduce the following random times tq = 0,(T„ = inf{t > T„,Xf € C},Tn = 
inf{t > cj„_i,Xf G 9(7}. We consider the Markov chain Z„ = X^ for n > 0. We see 
that from n = 1 on Z„ G dg. Also, X^^ can be any point of C, all the following X^^ 
belong to one of the Fj's. The chain never stops. 

We are now ready to prove: 

Theorem 3.2. For any 7 > there exists po > (which can be chosen arbitrary 
small) such that for any P2, < p2 < po, there exists pi, < pi < p2 such that for all 
X in the p2 -neighborhood of Oi{i = l,...,l) the one-step transition probabilities of Zn, 
Zq = X satisfy the inequality 

exp[-e-2(y+(0,, O,) + 7)] < P(x, dg,) < exp[-e-\V+{Oi, O,) - 7)] 

for some < e < Eq. 

Proof. We can assume V^{0i,0j) < 00. Set V^^ = max V^{0i,0j). Choose 

i,j=l,2,...,l 

positive Po small enough. For every pair Oi, Oj for which V^{Oi, Oj) < 00 we choose a 
function c/?^'-' G C [o,t] ( [-D] ) > < t < T = T(0i,0j), such that c/Jq-' = Oi,iflp = Oj,ipl'-' 
does not touch |J {Os}, and such that (by Lemma 3.2) 

S^rri^'^^) <V^{Oi,Oj) + 0.5^ . 
We choose positive pi smaller than ^, p2 and 

^mm{p{^l'\ \J{Os}):0<t<T,i,j = l,2,...l} . 

For every x in a /92-iieighborhood of Oj we take a curve connecting x with Oj and 
for which the value of S*"*" does not exceed 0.87 (by Lemma 3.1). We combine this curve 



with the curve c^J'-' and obtain a function ft,Ol^t<T,ipQ = x, ipT = Oj (with a possible 
small change of T from T = T{Oi, Oj)) such that 

S+r{^)<V+{O^,Oj) + 0.8-/ . 

Prom Lemma 3.2 we choose a Tq > T, and extend the curve ft to T < t < Tq hy 
using a trajectory of the dynamical system xt = b{xt) on diD, without changing the 
value of Sqj, (ip) from that of Sffj,{ip). We choose positive 6 less than pi,po — p2- For 
a trajectory of Xf starting from x, passing at a distance from ipt smaller than 5 for 
< t < Tq, it must intersect with Fj and reaches the 5- neighborhood of Oj without 
getting closer than 5 from any of the other Og, s ^ i,j. Moreover, X^_^ G dgj, thus 

P{x,dgj) > PApOToiX',^) <S}> exp[-e-2(5o+,,,((/.)+0.l7)] > exp[-e-^{V+{Oi,0,)+j)] 

Now we turn to the proof of the upper estimates. For any curve (pt,0 < t < T 
beginning at x, touching the ^-neighborhood of dgj, not touching any of the Og, s ^ i,j, 
we have 

S+riif) >V+{0^,Oj) -0.7 J . 

We use Theorem 3.1 to choose Ti such that for any x € [D] \ g we have P|{ti > 
Ti} < exp(-e-2y+) for some V^^ > 0. 

Any trajectory X^ beginning at x and being in dgj at time ti either spends time 
Ti without touching dg or reaches dgj over time Ti, in this case 

poTAx',<^Av+{o^,Oj) - 0.77)) > 5 . 

Therefore we have 



< exp(-e-2^+) + exp[-e-2(y+(0„0,) - O.97)] 
<exp[-e-2(y+(0„0,)-7)] 

for sufficiently small e. D 

In an exactly similar way one can also formulate the estimate on transition proba- 
bility based on the quantities 



y+(x, O,) = inf {S+^{f) :fo = x,fT = Oj,ft e[D]\\\ {O,}, < t < T} 

</peC[o,T] ([£>]) ^. 

We have 



Theorem 3.3. For any 7 > there exists po > (which can he chosen arbitrary 
small) such that for any P2, < P2 < Po, there exists pi, < pi < p2 such that for all x 
outside the p2 -neighborhood of Oi{i = 1, ...,/) the one-step transition probabilities of Zn, 
Zq = X satisfy the inequality 

expl-e^^V+ix, Oj) + 7)] < P(x, dg,) < exp[-e-\V+{x, O,) - 7)] 

for some < e < Eq. 

3.3 The invariant measure of X^; sublimiting distribution 

In this section we study the invariant measure of the process X^. Based on the 
estimates on transition probabihties given above, the proof of the asymptotic result is 
the same as that of [Sj Ch.6] and [1]. Let us formulate and prove two more technical 
lemmas, after which the rest of the proof is just a study of Markov chains on graphs. 
The latter part will be omitted since it is the same as O Ch.6] and [3]. 

Lemma 3.3. Fori G {1,2,...,/}, define 

T.^^O,) = mf{t,Xo' = x,Xf G dSsiO,)} . 

For any 7 > 0, there exist 6 > such that for all sufficiently small e and x € SsiOi) we 
have 

K-^Ss(o.) < exp(7e"^) • 

Proof. Choose point z £ D close to Oj. Put S = 2 • Connect x with Oi and 
Oi with z with the values of S^ not exceeding j and ^, the resulting function is called 
iff The length of the time interval of ipt is uniformly bounded by Tq for all x £ G. We 
extend ipt up to Tq by using a trajectory of xt = b{xt) in D U dD without making S^ 
larger. 

Now we have for x G SsiOi), 

Kire.m < To} > Pl{poToiX',ip) < 5} > exp(-0.97e-2) . 
Using the Markov property we see that 

K{r£,m > nn} < [1 - exp(-0.97£-')]" . 

This yields 

00 
Krssm < nJ2[l - exp(-0.97e-')]" = Tq exp(0.97e-') . 



n=0 



10 



Sacrificing O.I7 in order to get rid of Tq we get the desired result. D 

Lemma 3.4. For any 7 > there exist pi > such that for all sufficiently small 
£ and y G dgi we have 



/ XgA^!)dt > exp{--f£- 
Jo 



y , 
/o 



Proof. Choose pi smah. We connect y G dgi with Oi using a curve (ft, extend 
it using the trajectory of xt = h{xt) on dD till first exit time ctq from £pg{Oi), with 
corresponding S~^ less than O.57. All the trajectories at a distance less than ^ spends 
a time at least to > within gi, uniformly for all y € dgi. The probability of all such 
trajectories is no less than exp(— 0.97e~^). Thus the expected value is no less than 
to exp(— 0.976^^). By sacrificing O.I7 we can get rid of to- □ 

The rest of this section is devoted to the description of the algorithm for the calcu- 
lation of the invariant measure and the metastable states. The proof we shall omit here 
follows [a Ch.6] and [1]. 

Let L be a finite set (in our case L = {1,2, ...,l}), whose elements are denoted 
by letters i,j,k,m,n, etc. Let a subset W be selected in L. A graph consisting of 
arrows m —^ n (m G L\W, n G L,n ^ m) is called a M^-graph if it satisfies the following 
conditions: 

(1) every point m € L\W is the initial point of exactly one arrow; 

(2) there are no cycles in the graph. 

Intuitively, a M^-graph is a graph consisting of arrows starting from each point 
m G L\W , and going along a sequence of arrows leading to some point n G W . 

The set of VF-graphs is denoted by GiW). We shall use the letter g to denote 
graphs. 

Let W{Oi) = min Yl{m^n)eg'^^ iOm,On)- It can be proved that 

W{Oi)= min V V+{0m,0n) ■ 



9GG{i} ^ 



We have 



e 



Theorem 3.4. Let p,^ he the normalized invariant measure of the process X 
Then for any 7 > there exists pi > such that we have 

exp[-e-2(T^(Oi) -miniy(Oi) +7)] < M"(5i) < exp[-e-2(Vr(Oi) - min P7(0i) - 7)] 

i i 

for sufficiently small e > 0. 

11 



We shall say that a set N C [D] is stable if for any x £ N, y ^ N we have 
V~^{x, y) > 0. One can show that for an unstable Oj (j = 1, ..., /) there exist a stable Oj 
(i/ j,i = 1,...,/) such that y+(Oi,Oj) = 0. 

Theorem 3.5. For x G [D] set 

W{x) = mm[W{Oi) + V{Oi,x)] , 

where the minimum can be taken over either all of Oi, ..., O; or only stable ones. Let fi^ 
be the normalized invariant measure of the process Xf . Then for any 7 > there exists 
p > such that for any < p < p we have 

e^p[-e-'^{W{x) - mmW{Oi) + j)] < n%£p{x)) < eicp[-e-'^{W{x) - m.mW{Oi) - -f)] 

i i 

for sufficiently small e > 0. 

Here Sp{x) is a /j-neighborhood of x. 

The above two theorems roughly say that as first t ^ 00 and then e ^ 0, the 
process Xf will be situated in one of the Oj's which minimizes the values of W{Oi) (it 
can be calculated either via all Oi,...,Oi or only via the stable ones). In generic case, 
when min W{Oi) is attained at some unique point i, we have for any 6 > 0, 

i 

limlimP^{|Xf-O,|>5}=0. (3.1) 

e—^Ot—>-co 

A natural question is that how the limiting distribution behaves when we take 
the limit in a coordinated way, i.e. take e — )• and t = t(e^^) — > 00. This is the 
problem of metastability and sublimiting distributions (see [2]). Let us assume that 
T = T{e) X exp(^) and we consider limP|{X|,.^-, G F}. In the generic case one can 
define a function K*{x, A) € {1, 2, ..., Z} such that 



limP^{|Xf(,) - 0;^,(,,A)l > n = (3.2) 

for any 5 > 0. 

The algorithm to determine K*{x, A) is as follows. First we consider for each Oi (the 
rank cycle) the "next" most probable w-limit set J\f{Oi) that we are going to jump to. 
Continuing this determination of "next" states we form the rank 1 cycle Oj -^ J\f{Oi) -^ 
M^{Oi) -^ ... -^ ATi^iiOi). We stop once we get a repetition AA(7\r'0}(Oi)) = Oi. 
Cycles generated by distinct initial points i £ {!,...,/} either do not intersect each other 
or coincide: in the latter case the cycle order on them is one and the same. 
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We continue by recurrence. Let the cycles of rank (A;— 1) be vr-^^" ,...,7r^r^ . Starting 
from each (/c — l)-cycle vr^ ~ one can determine the "next" most probable [k — l)-cycle 
TT ~ i. 1, that we will first iump to. Continuing this determination we form a rank 

k cycle tTj ~ — )■ it 7 ^_j^ —> ... ^ ii ~m ,,_^ . We stop once we get a repetition 

M{N "i {tt- )) = T^i ■ Cycles of rank k generated by distinct cycles of rank k — 1 
either do not intersect each other or coincide. 

In this way we can continue until the last cycle which is the whole of {Oi, ..., O;}. 
The metastable states are determined by the timescale of the cycles that we traverse. 

Let us be more precise. Starting from a cycle vr, to determine the "next" cycle 
M{tt) that we first jump to, we calculate 

A{7t)= min V V+{Ora,On) . (3.3) 

geG{L\7r) ^ 

Here L = {1,2, ...,/}. The minimum of the above expression determines a L\k 
graph consisting of chains of arrows leading to the first state in L\7r we jump to. 
We put 

C(7r) = A(7r) - min min V F+(0„„0„). (3.4) 

[m^rnjag 

Here G-,^{i} is the set of {i}-graphs restricted to vr. Then the asymptotic exit time from 
vr is of order x exp ( — ^ J . 

Starting from i = i{x) (which is the label for the first equilibrium among Oi, ..., Oi 
that we approach in finite time, starting from x), let vr, vr',..., vr^''-' be cycles of next 
to the last rank, unified into the last cycle, which exhausts {1,2, ...,/}. If the con- 
stant A is greater than C(7r), C(7r'), ..., C(7r*-''-'), then over time of order exp(Ae^^) the 
process can traverse all these cycles many times (and all cycles of smaller rank inside 
them) and the limiting distribution is concentrated on that one of the cycles for which 
C{'k), C(7r'), ..., C{'K^'^') is the greatest. Within this cycle, it is concentrated on that one 
of the subcycles for which the corresponding constant C(») in (3.4) is the greatest possi- 
ble, and so on up to points (one point in the generic case) Ox*(x,A)- This point Ok*{x,\) 
is the metastable state in (3.2). 

4 Application to PDE 

The solution of (1.1) can be represented through process (1.2) by the formula 
u^(x,t) = 'Eixg{Xf). This is an immediate consequence of the following generalized Ito's 
formula: 
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Lemma 4.1. Assume process {Xf,(^f) is given by (1-2), Xq = x. Let u{x,t) be of 
class C^'^(M X R_|_) with uniform bounded derivatives up to the second order in x and 
up to the first order in t. Then we have 

u{XI,t)-u{x,0) 

= ^ (^^ + C'^ u{Xl s)ds + ^ Vu{Xl s) ■ 7(Xf )C + j^ Vn(X|, s) • a{Xl)dWs . 

For a proof of this theorem see ^ Section 3] . 

Our answer to the problem (1.1) is 

Theorem 4.1. Under all our assumptions, in generic case, for T{e) x exp(T), 
we have 



linin^(x,T(e)) = g{OK*{x,X)) 



where K*{x,X) is defined as in Section 3.3. 



5 Example 



Consider an example. Let the domain Z) be a unit disk -B(l) = {{yi,y2)',yi + 
^2 < 1} in I^^. Let the smooth vector field by{yi,y2) be given such that by{yi,y2) = 
(byi{yi,y2),by2{yi,y2)) is as in Fig.l. We consider the problem 



du''{yi,y2,t) e 



2 



■^yi,y2u''{yi,y2,t) + by{yi,y2) ■Vu^{yi,y2,t) , s>0 



dt 2 

u^{yi,y2,0) = g{yi,y2) , y? + yi<i; 

— {yi,y2,t) = , yf + y^ = l,t>0. 

(5.1) 

Here 7— is the derivative with respect to the inward unit normal. The action 
or 
functional takes the form 



' 1 /"^ 

+ I - / \(ps - by{(ps)\K2ds , (f G C[o,T]([-D]) absolutely continuous ,(fo = x ; 

'-'oTvV'i ~ \ -^ Jo 

+00 , for the rest of 99 G C[ot]([-^]) • 

(5.2) 

We calculate the "quasi-potential" using (5.2) 
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V+{x,y)= inf {S^rj.{Lp),(po = x,ipT = y, (ft & D U dD,0 < t <T < 00} . 

VGC[i,,T]([0]) 

The (jj-liniit sets of the dynamical system it = hy{xt) are the zeros of the vector 
field hy{x) on dD = S^. (And also the origin but it is unstable so that we neglect it.) 
In Fig.l the points Oi, O3 and O5 are stable equilibriums are the points 02, O4 and Oq 
are unstable ones. We can consider only the quasi-potentials between the stable ones. 
Suppose we have ^+(01,03) = 1, V+{03,0i) = 2, V+iOuO^) = 6, V+{05,0i) = 7, 
V+{05,Os) = 3,V+{03,05)=4. 

We are concerned with the limit lim.u^{yi,y2,T{e)) for T{e) x exp(^). Starting 

£4,0 ^ 

from the initial point (yi, y2), we suppose that we are attracted to Oi first. By calculating 

min Yl V^{Om, On) = 1 we see that over time exp(p-) we are going to jump 

geG{L\{l}) („^„)gjj 

to O3 first. We then calculate min ^ V^{0m,0n) = 2 and we see that over 

geG(L\{3}) (^rn^n)£g 

time exp(^) we will jump from O3 back to Oi and we form a cycle vr^^' = {1,3} of 

rank 1. We then calculate ^(7r(i)) = min ^ y+(0™, 0„) = F+(Oi, O3) + 

<yeG(L\7r(i))(„j^„)gg 

^+(03,05) = 5 and the first state out of cycle vr'-*^' that we are going to jump to 

is O5. Within cycle ir^"^' we are mostly staying in O3. We calculate C{tt^'^') = 5 — 

min min ^ y+(Om,0„) = 4. This means, that over time exp(^) we are 

jumping from O3 to O5. We then calculate min ^ y+(Om,0„) = 3 and we 

3eG(L\{5}) („„^„)g^, 

see that we are jumping from O5 out to O3 in time exp(-|-). This implies that within 
the cycle vr'^-' = {1,3,5} which exhausts all tj-limit sets, we are mostly staying in vr'-"'^', 
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and within t:^^' it is O3. 

Our result can be summarized as 

limu'^(yi,y2,T(e)) = giOi) for T{e) x exp(^) and < A < 1 ; 

l[mu''{yi,y2,T{e)) = giOs) for T(e) x exp(^) and 1 < A . 

Acknowledgement: We would like to thank our advisor Professor Mark Freidlin 
for posing this problem to us and for many useful discussions. 
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